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Abstract. We show how to find series expansions for vr of tlie form 
vr = X^Ji^o '^('^)/(^n)'^"' where S{n) is some polynomial in n (depending on 
m,p,a). We prove that there exist such expansions for m = 8k, p = Ak, 
a = (—4)'^, for any k, and give explicit examples for such expansions for small 
values of m, p and o. 



1. Introduction 
Using the formula (due to Bill Gosper 0) 

°° 50n — 6 



TT 



\n)^ 



Fabrice Bellard file pil.c] found an algorithm for computing the n-th decimal of 
vr without computing the earlier ones. Thus he improved an earlier algorithm due to 
Simon Plouffe 0. 

This formula can be proved in the following way. We have 



Hence the right hand side of the formula will be 



xTdx. 



t-l °° / 2(-\ ^X" 

/ V(50n-6)(3n+l) ( "^ U " 
-^0 n=o V 2 y 
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But 



so we get 



^(50n-6)(3n + l)y = — -— , 

n=0 ^ 



Z.OX — uux -r ^ox — y/x -r »/x — u 
RHS — o — „ ax 

Jo 



28x^ - 56x^ + 28x^ - 97x^ + 97x'^ - 6 

- x2 + 2)3 

4:x{x - l){x^ - 2Sx^ + 9x + S) 



+ 4 arctanfa; — 1] 



= TT. 





(^3 _ ^2 + 2)2 

We then asked if there are other such formulas of type 

S{n) 
2^ Ma-' 

n=0 \pnJ 

where S{n) is a polynomial in n. Using the same trick as before we have 

-i^ = (mn + l) f x^-{l-xf--^>dx, 

\vn) Jo 



and hence 



RHS = J ^(mn + l)S{n) ^ ^' j dx. 

If S{n) has degree d then 



(1.1) 



n=0 ^ ^' 

where T has degree d-\-l. It follows with 

xP(l - xY'^P 

y = 

a 

that ^ 

i?i/5= /" -— ^ UT^^dx, (1.3) 

where P{x) is a polynomial in x of degree m(d+ 1). We want this integral to be equal 
to vr. A good way to get tt is to have arctan(x) or arctan(x — 1) after integration. 
This means that 

x^\\- x)"'-^ - a 
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must have the factor + 1 or (x — 1)^ + 1, that is have a zero at i or This restricts 
m and p and gives the value of a. After experimenting with the LLL-algorithm we 
found formulas for vr in the following cases (there could be many more): 



m 


P 


a 


deg(S) 


3 


1 


2 


1 


7 


2 


2 


5 


8 


4 


-4 


4 


10 


4 


4 


8 


12 


4 


-4 


8 


16 


8 


16 


8 


24 


12 


-64 


12 


32 


16 


256 


16 


40 


20 


_45 


20 


48 


24 


46 


24 


56 


28 


-47 


28 


64 


32 


48 


32 


72 


36 


-49 


36 


80 


40 


410 


40 



Gosper 0| 
Bellard 



E.g., to find the formula in the last case we computed 

00 ^ 

n=0 V40n/ 

for k = 0,1,2, ...,40 with 6000 digits. Then the LLL-algorithm found (in about 
three days) the linear combination between vr and s(0), s(l), . . . , s(40). A good check 
is that the coefficient of vr is a product of small primes (actually all less than 80). 



2. Proving the formulas 



For m ^ 16 the integrals can be computed by brute force using Maple. The higher 
cases are all symmetric, i.e., m = 2p. Using the symmetry 

X < — > 1 — x 

^ This formula was used in Bollard's world record setting computation of the 1000 billionth 
binary digit of tt, being based on the algorithm in [Q. 
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or 



one can assume that arctan(x) and arctan(x — 1) each contribute 7r/2. So we make 
the wild assumption that 

P(x) , R(x) , , , , 

, ax = , + 2 arctanfa;) + 2 arctanfx — 1), 
Q[x)"-~^'^ Q[x)'^^^ 

where R{x) is a polynomial with R{0) = -R(l) = 0. Differentiation with respect to x 
yields 

P _ R! (A ]]'^'^ of ^ ^ 

QR' - a + DQ'R = P - 2Q-^ + ^dlTl) ■ 

This is a differential equation of first order in R with polynomial coefficients. Taking 
the case ('^g"), we have deg(P) = 544 and d = 16, so the denominator Q'^^'^ has 
degree 576. We solve the differential equation using Maple's 

dsolve(*,R(0)=0, series) 

setting the order to 600. Then we find that i? is a polynomial of degree 543. In 
practice one replaces the coefficients of -^^^^ and ^i_l^_^_2 by 2 times the denominator 
in the formula found for vr. Then R will have (huge) integer coefficients. 

Remark 1 . Computing the integral in the case (^^) , Maple got the terms 

/2x3 - + 7x - 3\ I 
arctan + arctan 

V 5 ; V2 4 

which can be shown to be equal to 

arctan(x) + arctan(x — 1). 

In the nonsymmetric case {^^^) one gets arctan(x — 1) and arctan(x'^ — 2x^ + x — 1), 
but the latter term makes no contribution since it has the same values at and 1. 

Remark 2. When we set this up then, in principle, we could let Maple do the eval- 
uation of T{y) (and thus of P{x)) through ( p.. 21) , i.e., let Maple's formal summation 
routines compute the sum 



n=0 
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or, to make the point clearly visible, compute a sum of the form 

oo 



n=0 



where S{n) is a certain polynomial in n. However, it turns out that this is very 
time-consuming. It is a much better idea to expand S{n) as S{n) = J2'j=o 

and then use that J2n=o G+D^" = - vY^^ for j = 0, 1, . . . , + 1. See also 

the approach that we follow in Section when we prove our main result. 



3. The CASE (^^^ ), = 1,2,3, ... 

Solving the differential equation for R{x) one finds that R{x) has the factor 
a;(l — x){2x — 1). If one divides out the factor 2x — 1, the remaining factor is invariant 
under the substitution x < — ^ 1 — x. Hence we make the substitution 



where 
Then 

Now 

so 

and 



R{x) = (2x- 
z = x{l — x). 



dz 

— = 1 - 2a; and (1 - 2xy = 1 - 4z. 
ax 



Q = z^'- (-4)^ 

Q' = Akz^''-\l - 2x) 



R! = 2R{z) + {2x - 1)^(1 - 2x) = 2R - {1 - 4z)^. 

dz dz 



It follows that the left hand side of the differential equation is 

QR' - (4fc + 1)Q'R = q(2R-{1- 4^)^^ + (4A; + l)4A;z^'=-^(l - 2xfR 



dR 

(2Q + 4k{4k + 1)(1 - 4z)z^''-^)R - (1 - 4z)Q^. 

dz 
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On the right hand side we have 



1 1 3-22 

+ 



and Q is for all k divisible by 

+ A= {z'^ -2z + 2){z^ + 2z + 2). 

Hence P and Q and 



are all polynomials in z. Finally we get the following differential equation for R: 

-{l-4z)Q — + {2Q+4k{4k+l){l-4z)z"'-')R-P+2{3-2z)^^ = 0. (3.1) 

If P is known then one solves this equation for R just as before (but the degrees 
are cut in half). 

But there is a possibility to find the formula for vr and 'prove it in one stroke. Let 



and assume that 



and 



N{k) = 4A;(4fc + 1), 



N{k)-l 

R{z)= J2 «0>' (3-2) 



Skjn) , . 

n=0 \Akn) V ^1 



with 



4fc 

S,{n) = Y,<N{k)+j)n^. (3.4) 
i=o 

Substituting this into the differential equation we will get a system of linear equations 
for the a(j)'s of size {N{k) +4k) x {N{k) +4fc). If this system is nonsingular, we can 
be sure to be able to solve it and thus find and prove an expansion for tt of the form 
(P). 

We went (again) to the computer and generated the system of equations for small 
values of k. Let A{k) be the matrix of the system. Let further r{k) be the least 
common denominator of the coefficients of Sk{n). We obtained the following tables: 
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1 
2 
3 
4 
5 



det{A{k)) 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 



291385772 

-2523352517714114133 
21367317754172511201319175194232 
.232313I67583753112813271725198236293312 

2539932905345793ii41;L337l733l93223l02973l6373 

r(fc) 



325272 

365372112132 

2^335272112132172192232 
23310567311 . 132172192232292312 

2113135374113133172192232292312372 

210375873113133172192232292312372412432472 

2173165973114133173192232292312372412432472532 

2103195878114133173193232292312372412432472532592512 

2233451078114134173193233292312372412432472532592512572712 

2203205479115134173193232292312372412432472532592512572712732792 



From the table it seems "obvious" that the determinant of A{k) will never vanish 
as it grows rather quickly in absolute value with k. We also see that the same primes 
occur in the factorizations of r{k) and det{A{k)). Even more striking, it seems that 
the largest prime factors occurring grow only slowly (namely approximately linearly) 
when k increases. The last fact strongly indicates that there may even be a closed 
form formula for det{A{k)). As it turns out, this is indeed the case. In the next section 
we will explicitly compute the determinant of an equivalent system of equations, from 
which it follows that (see the remark after Theorem |]) 



=((4^ + i)')-M!nM 



det{A{k)) = {-iy-'2 



fc-lo32A:3+24A;2+2A:-1^8fc2+2fc 



(3.5) 



Hence we have the following theorem, which is the main theorem of our paper. 
Theorem 1 . For all k > 1 there is a formula 

Sk{n) 



; (a(-4) 



kn ' 



where Sk{n) is a polynomial in n of degree 4k with rational coefEcients. The polyno 
mial Sk (n) can be found by solving the system of linear equations generated by (PTT 
and the Ansatz (|3.2| ) and (^l4| 
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The denominators of all a(j)'s divide r{k) which is much smaller than det{A{k)). 
This means that there must be some miracle occurring at the end when solving the 
system. E.g., when k = 5 then det(A(5)) has about 2400 digits but r(5) only 40. 
Unfortunately, we are not able to offer an explanation for that. 

In practice we are only interested in the coefficients of Sk so we try to eliminate 
a{l) , a{2) , . . . , a{N (k) — 1) first. This can be done by first avoiding all equations 
containing a{N{k)), . . . , a{N{k) +4A;), i.e., the equations coming from the coefficients 
of t^'^'", V = 0,1, ... ,4k + 1. This looks very nice theoretically, in particular as the 
system for a(l), . . . , a{N{k)) is triangular, but in practice the computer breaks down 
since the rational numbers occuring become very large. 

We close this section by listing a few explicit examples. 
Example 2. We have 

TT 

where 
and 

S{n) = -89286 + 3875948n - 3497013471^ + 110202472n=^ - 115193600n^ 
Example 3. We have 




where 

r = 3^5^7211^13^ 

and 

S{n) = -869897157255 - 3524219363487888n + 1124667772631 18189/1^ 
- 1242789726208374386n^ + 6693196178751930680n^ - 19768094496651298112n^ 
+32808347163463348736n^-28892659596072587264n^+10530503748472012800n^ 



S{n) 



r = 3^527^ 



Example 4. We have 

_ 1 ^ Sin) 



^ n (?«")256"' 

n=0 \16n/ 

where 

r = 2^3^°5^7^11 ■ 13217219^232292312 
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S{n) = - 2062111884756347479085709280875 

+ 1505491740302839023753569717261882091900n 

- 11240114940408765821383938671621 1975291975n2 
+ 3257881651942682891818557726225840674110002n^ 

- 51677309510890630500607898599463036267961280n^ 
+ 517337977987354819322786909541179043148522720n^ 
-3526396494329560718758086392841258152390245120n^ 
+ 171145766235995166227501216110074805943799363584n^ 
-60739416613228219940886539658145904402068029440n^ 
+ 159935882563435860391195903248596461569183580160n'^ 
-313951952615028230229958218839819183812205608960n^' 
+ 457341091673257198565533286493831205566468325376n^ 
-486846784774707448105420279985074159657397780480n^' 
+ 367314505118245777241612044490633887668208926720n^" 

- 185647326591648164598342857319777582801297080320n^' 
+ 56224688035707015687999128994324690418467340288n^^ 
-7687255778816557786073977795149360408612044800n^^ 



Example 5. (nonsynunetric). We have 




where 

r = 3^5 • 7^11 • 13 • 17 • 19 

and 

S{n) = -4843934523072 - 1008341 177146848n + 23756198610834352n^ 

- 242873913552020704n^ + 1195813551184400032n^ - 3272960363556054592n^ 
+ 490937916783701 1328n^ - 3816399750842818816n^ + 1190182007407360000n^. 
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4. Proof of the Theorem 



We want to prove that, by making the Ansatz ( |3.2| ) and ( p.4| ) and substituting this 
into the differential equation ( |3.1j ) (the polynomial P being given by Sk through ( |1.1|) - 



(|1.3|), Q being given by (x(l — x)Y'^ — (—4)'^), the resulting system of linear equations 
will always have a solution. In fact, we aim at finding an explicit formula for the 
determinant of the corresponding matrix of coefficients that allows us to conclude 
that it can never vanish. 

It turns out that for that purpose it is more convenient to set up the system 
of linear equations in a different, but equivalent way. This equivalent system will 
have coefficient matrix M (see ( [4.4|) ). The evaluation of its determinant will be 
accomplished through Eqs. (|0| ), (g3)> (ED. (EH), and Theorem |. 

To be precise, we encode the polynomial Sk{x) (and, thus, T{y) and P{x)) differ- 
ently. We claim that T{y) has an expansion of the form 

ny) = c,y\ (4.1) 

subject to the single constraint 

Y[(4zk-l)co+Y.{-iy I n m-l)]ll[{Mk + l)]c,=0. (4.2) 
i=i j=i \ i=i J \i=i J 

(As usual, empty poducts have to be interpreted as 1.) This is seen as follows. The 
polynomial Sk{x) can be written in the form 



j=0 ^ / 
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for some coefficients sj. Hence, we have 
-^M_ = f;(4A:n + l)5.W 

^ n=0 

n=0 jr=0 ^ ^ 

= E + J + 1) - 4A:(J + 1) + 1) (''4^ ^" 

j=0 n=0 ^ ^ 

4fc 00 / / _|_ ■ _|_ 1 \ 



j=0 n=0 
Ak 



4:k r 



= E + 1) (Y^^ - mJ + 1) - 1) (T^^ j 

^ 4A:+1 

= n- )4fe+2 E + 1)^4^-, + (4A;(4fc - J + 2) - l).4.-i+i) • 

In the last hne, s_i and S4A:+i have to be read as 0. It is now a trivial exercise to 
substitute the coefficients of in the sum in the last line into the left-hand side of 
and verify the truth of ( |4.2|) . 
The above implies that 

Ak+l 
j=0 



where the coefficients Cj obey 

Now we are ready to set up the system of linear equations. We make again the 
Ansatz (|3.2|) , but we replace ( |3.4| ) by 

Ak+l 

P{z) = <N{k) + j)z'^'/{-^y', (4.3) 

j=0 

where the a{N{k) +j), j = 0,1, . . . ,4k + l, are subject to (i.e., the relation 
holds when Cj is replaced by a{N{k) + j)). Clearly, we have to add (^?^ ) to the set of 
equations that result from the differential equation (3.1). 
The coefficient matrix of the system looks a follows: 
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where x is a line vector of N{k) — 1 zeroes, y = {yo,yi, ■ ■ ■ ,1/4/c+i) is the vector of 
coefficients of ( [4 .21) , i.e., yo = Yltti^i^'^^ " l)' 



£ = 1, 2, . . . , 4fc + 1, where U is an {N{k) + ik) x {N{k) - 1) matrix and V is an 
{N{k) + 4k) X [Ak + 2) matrix, both of which we define below. 

We consider the top-most line of M (which is formed out of x and y) as row of 
M. We label the rows of U and V hj i running from 1 to N{k) + 4k. Furthermore, 
we label the columns of M by j running from 1 to N{k) + 4k + 1. 

Following this labelling scheme, the matrix U has nonzero entries only in the four 
diagonals i=j,i = j + l, i = j + 4k, i = j + 4k + 1. We denote the entries in 
column j on these four diagonals in order /o(j), /i(j), go{j), and gi{j), where 



To be precise, the (j, j)-entry is /o(j), the (j + 1, j)-entry is the (j + 4A;, j)-entry 

is goU), the (j +4k + 1, j)-entry is gi{j), j = 1, 2, . . . , N{k) - 1. 



On the other hand, the matrix V is composed out of columns, labelled N{k), 
N{k) + 1, . . . , N{k) + 4A; + 1, each of which containing just one nonzero entry. To 
be precise, the nonzero entry of column N{k) + j is located in the {4jk + l)-st row 
(according to our labelling scheme), and it is equal to (— 4)"-'*^, j = 0,1, . . . ,4k + 1. 

We will now compute the determinant of M and show that it does not vanish for 
any k. 

We perform some row operations on M, with the effect that the entries of y get 
eliminated. This is achieved by subtracting ni=t^(4^^ ~ 1) times row 1 from row 0, 
and 



times row 4£k+l from row 0, £ = 1, 2, . . . , 4A; + 1. Doing this, we must expect changes 
in row in columns 1, 4k, 4k + 1, 8A;, 8A; + 1, ... ,N{k) —4k = 16/c^. However, at this 
point a miracle occurs: the new entries in row in columns 4ik+l, £ = 0, 1, ... , 4k — 1, 
are still 0. On the other hand, the values of the new entries in row in columns 4£k, 




goU) 
giU) 



J(-4)^ 

-(4J + 2)(-4)^ 

iN{k)-j), 

~{4N{k)-4j-2). 
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1,2, 



, Ak are 



-11 



(4.5) 



After these manipulations we obtain a matrix of the form 

M' -- 



X 

U 



y 

V 



where x' and y' have the same dimensions as before x and y, respectively, where the 
only nonzero entries of x' are in columns labelled by numbers which are divisible by 
4k, with the entry in column 4ik given by ( ^.5| ), and where y' consists only of zeroes. 
We have 

detM = detM'. (4.6) 

The next step consists in expanding the determinant of M' with respect to columns 
N{k), N{k), . . . , N{k) + 4fc + 1 (i.e., the last 4k + 2 columns). Since each of these 
columns contains just one nonzero entry (which is a power of (— 4)"^^), we have 

det M' = ±(_4)-fe(2fc+i)(4fc+i) ^//^ 

where M" is the matrix arising from M' by deleting the last 4k + 2 columns and the 
rows 4ik + 1, £ = 0,1,...,4A; + 1. More precisely, the matrix M" has the following 
form: 



/ Xi 


X2 




Fi 








Gi 


F2 








G2 


Fz 















V 



Gik 












F^k+i 



where x^, £ = 1, 2, . . . , 4k, is a line vector with 4k entries, all of them being zero except 
for the last, which is equal to (|4.5| ), where x^k+i is a line vector of 4/c — 1 zeroes, where 
and G^, £ = 1,2,..., 4k, are {4k — 1) x (4A;) matrices with nonzero entries only in 
the two main diagonals, and where -Fl^fc+i and G^k+i are {4k — 1) x {4k — 1) matrices, 
Gi^k+i being upper triangular. To be precise, for £ = 1, 2, . . . , 4A; we have 

//i{4{<?~l)fc+l) /o(4(^-l)fc+2) ... \ 

' /i(4(^-l)fc+2) /o(4(^-l)fe+3) 



/i(4£fc-2) /o(4£fc-l) 

h{Uk~l) h{Atk) / 
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and 



/fli(4(€-l)fc+l) go(4(f-l)fe+2) 
' fli(4(£-l)fc+2) so(4(£-l)fc+3) 



and we have 



ai(4£fc-2) go(4ffe-l) 

gi{4ek-l) go{ilk) J 



/fli(16fe2+l) ao(16fc2+2) 

fli(16A;2+2) go(16A;2+3) 



4fe+l 



V 



fli(16fc2+4fc-2) go(16fc2+4fe-l) 

gi(16fc2+4A;-l) / 



The precise form of F^k+i is without relevance for us. Wc do a Laplace expansion 
with respect to the last 4A; — 1 rows. Because of the triangular form of G^k+i we 
obtain 

'16fc2+4fc-l 



detM"= ( W gi{i) ) detM'", 

=16fc2+l 



(4.8) 



where 



M" 



( Xx X2 X3 

Fi 

Gi F2 

G2 F3 

G3 



V 





G4fe_i F4k 

G4fe/ 



Instead of M'" we consider a more general matrix. Define the functions 
Mt,j)^((N(k)+j)Yt-X2,t)(-^)\ 

Mt,j) = -amk) + 4j + 2)y, - 4Xl,,)(-4)^ 

goit,j)^{X2,t-jYt), 



It should be noted that these functions specialize to fo{j), fi{j), go{j), gi{j), respec- 
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tively, if Xi^t — -^^2,* = N{k) and Yt — 1. Now we define the matrix by 



M 



X 



X2 


X2. 





















F,^ ... 














r^x 












(4.9) 



where 



//i(l,4(^-l)fc+l) /o(l,4(^-l)fc+2) 







/l(2,4(^-l)fc+2) /o(2,4(^-l)fc+3) 



/i(4fe-2,4tt-2) /o(4ifc-2,4^fe-l) 

/i(4it-l,4^fe-l) /o(4ik-l,4^fe) / 



and 



/gi(l,4(<?-l)fc+l) so(l,4(^-l)fe+2) 

gi(2,4(^-l)A;+2) so(2,4(^-l)fe+3) 



v 



gi{ik-2,Aik-2) go(4fc-2,4£fe-l) 

gi(4fe-l,4tt-l) ffo(4A;-l,4^fe) / 



(4.10) 



Clearly, we have 

The evaluation of det is given in the theorem below. Prom the result it is 
obvious that det L ^ i is nonzero, and, thus, also det M. 



4k 



Theorem 6. We have 

det = (_l)A.-l42fe(4fe^+7fe+2)^2fe(4;.+l) -Q^^ ^ 

i=l 

4A:-1 

X J] (2Xi,„ - (32A;2 + 2a - l)r„) 

a=l 

X J] (2X2,^1; -2Xi,„n- (26- 2a + l)F„n), (4.11) 



l<a<6<4fc-l 
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where {a)k is the standard notation for shifted factorials, {a)k '■= a(a + 1) ■ • • (a + 
k — I), k > 1, and (a)o := 1. 

Remark 3. Once having found this theorem, it is not difficuh to prove ( p.5|) , by 
working out how the coefficients of P{z) resulting from the Ansatz ( |4.1| )-( |0| ) are 
related to the coefficients of P{z) resulting from the Ansatz ( p.4|) . Since this is not 
essential for the proof of Theorem |T], we leave the details to the reader. 

Proof of Theorem ^. We follow the "identification of factors" method as described 
in Section 2.4 in 0. 

First we show that (2Xi,„ - {32k'^ + 2a - l)^^) divides det M^, a = 1, 2, . . . , 4fc- 
1. What has to be proved is that det M-^ vanishes for Xi^a = (32A;2 + 2a - l)Ya/2. 
This can be done by showing that for this choice of Xi ^ there is a nontrivial linear 
combination of the rows of . Indeed, if Xi^a = (32/c^ + 2a — 1)1^0/2 we have 

2(X2,4fc-l - (Nik) - l)y4fc-l) . ^ 

■ (row of M ) 



.4)fc(4fc+i)+i(i6A;2 + 1) nt=r(4tA; + 1) 

2^ I Ark 11 



4rfc lll6p+l_4tA; 

r=0 s=0 \ t=0 

^jT 2Xi,t - (32fc2 + 2t - l)yt 



• (row (16A;2 - (4A; - l)r - s - 1) of M^) = 0, 

as is easy to verify. 

Next we claim that {2X2,bYa-2Xi^aYb-{2b-2a+l)YaYb) divides det M^, 1 < a < 
6 < 4fc — 1. Let us ffist impose the additional restriction that a < b. Using the above 
reasoning, the claim then follows from the fact that if X2fi = T^Xi^a+ (26 — 2a + l)Yfc/2 



2Xi^a - (32A;2 - 8kt + 2a + l)Ya 
\ 2Xi^a - (64P + 8A; - Skt + 2a + 1)1^ 

2Xi,,y; - 2Xi,.y, - (2t - 2a)y.y, 

2X2,t^iYa - 2X^^aYt-i - {2t -2a- 1)^-1 




• (row (16A;^ - {4k - l)r - s - 1) of M^) = 0, 

as is again easy to verify. On the other hand, if a = 6, then the same argument 
shows that {2X2^a — '^Xi,a — Ya) divides detM"^. It remains to be checked that also 
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Ya divides det M^. Indeed, if Ya = then we have 

(row {IQk^ - {4k - l)r ~ 4k + a) of M 



Ak ^ 



r=0 ^ ^ 

whence Ya divides det for a = 1,2, ... ,4k — 1. 

These arguments show that the product on the right-hand side of ( |4.11| ) divides 
det polynomial in the Xi^^'s, X2,aS, and Y^'s. 

Clearly, the degree in the Xi^^'s, X2,aS, and YaS of det is at most 16/c^ — 1. 
But the degree of the right-hand side of ( 4.11| ) is exactly 16fc^ — 1. Therefore we have 
proved that 

4k 4fe-l 



det = C, Hit + l)4,.-i+i n (2^1,- - + 2a- 1)F,) 

X II {2X2,tYa-2Xi,aYt-{2b-2a + l)YaYt), (4.12) 



l<a<fe<4fe-l 



where Ci is a constant independent of the Xi ^'s, X2^aS, and F^'s. 
In order to determine Ci, we compare coefficients of 

X^Xj^r^ • ■ • Xl,_,Y,%' . . . Y,tzl (4.13) 



on both sides of ( 4.12 ). We claim that the coefficient of this monomial in detM"^ is 



equal to detM*-^, where M*-^ is defined exactly in the same way as (see ([4.9|) ), 
except that the definitions of the functions /o, fi,go,gi are replaced by 

fo{t,j) = {N{k)+j){-4)\ (4.14a) 

/l(^,J)=4(-4)^ (4.14b) 

9o{t,j) = -j, (4.14c) 

gi{t,j) = -4. (4.14d) 



This is seen as follows. The monomial ( [4.13|) does not contain any X2^a- Therefore 



for finding its coefficient in detM"^, we may set X2^a = in for all a. 

In which way may the monomial ( [4.13| ) appear in det (with all X2^a equal 



to 0)? A typical term in the expansion of det M'''- is the product of 16k^ entries 
of , each from a different row and column. The monomial ( 4.13 ) contains Xf\. 



The variable Xi^i is only found in columns 4ik + 1, i = 0,1, . . . ,4k — 1 (and rows 
labelled by numbers = 1 mod 4fc — 1, according to our labelling scheme). Therefore 
in a product of entries (each from a different row and column) which produces a 
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term containing Xf\ all the entries from columns 4ik + 1 must be ones containing 
Xi_i. This explains the above definitions ( [4.14b|) and ( [4.14d|) of /i(l,4£/c+ 1) and 



gi{l, 4:ik + 1), i = 0, 1, . . . , 4A; — 1, respectively. Moreover, we must generate the Yi 
in (|4.13|) from an entry in a column Aik + 2, for some i. (The variable Yi is also 
found in entries in columns Alk + 1, but these columns are already taken by our 
choice of entries which contain the Xi^i's.) This explains the definitions (|4.14a|) and 
(|4l4^ ) of /o(l, Aik + 2) and c/o(l, 4:ik + 2), £ = 0, 1, . . . , 4fc - 1, respectively. Next 



we ask how we can find (in the remaining columns and rows) entries which contain 
X^2~^- Arguing in an analogous manner, the variable X12 only appears in columns 
Aik + 2, £ = 0, 1, 2 . . . , 4fc — 1. One of these columns is already taken by the entry 
from which we picked Yi. Therefore in all the remaining ones we must choose entries 
containing Xi^2- This explains the definitions ( |4.14b| ) and ( [4.14dD of f i{2, Aik + 2) 



and gi{2, Aik + 2), i = 0,1, ... ,Ak — 1, respectively. Next we consider the term Y2 
in ( ^.131 ). It must come from two entries in columns Aik + 3, for two different €s. 
This explains the definitions ( ^A4a| ) and (|4A4cD of foi2,Aik + 3) and go{2,Aik + 3), 



i = 0,1, . . . ,Ak — 1, respectively. Etc. 

The evaluation of det follows from Lemma |^ below with Xa = 1 and Za = 
N{k) hi a = l,2,...,Ak-l. □ 

We consider now a more general determinant than det M'-' , the latter having been 
defined through the functions in ( |4.14| ). Replace these functions by 



fo{t,j) = {Zt+j){-A)\ 

f,{t,j)=A{-A)'X„ 

9o{t,j) = -j, 
gi{t,j) = -AXt. 

Let us denote the matrix defined by these functions in the same way as before M'-^ 
by . Clearly, specializes to M'" if all Xt are set equal to 1 and all Zt to N{k). 
The determinant of evaluates as follows. 

Lemma 7. We have 

4k-l / a-1 \ 

det = (_l)fc-l2l6fc^+20fc2+14A:-1^4A:^4^ + 1)' H ^a'^'"" H^^^ ' 4^^) " (^-IS) 

a=l \ 6=0 / 

Proof. We proceed in a similar way as in the proof of Theorem |^. In the first 
step we show that the product on the right-hand side of ( [4.15|) divides det as a 



polynomial in the X^'s and ZaS. Then, in the second step, we compare the degrees 
of the product and det . Since the degree of det turns out to be at most the 
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degree of the product, it then follows that det is equal to the product times some 
constant which is independent of the X^'s and ZfS. Finally, in the third step, this 
constant is found by computing the leading coefficient of detM^. 

Step 1. The product (Xf +1"'^ Hb^dl^a - 4&A;)) divides detM^. We start 

by applying several row and column operations to detM^, with the final goal of 
reducing the size of the determinant. First, for i — 16A;^ — 1, 16/c^ — 2, . . . ,4/c, in 
this order, we add (—4)*^ times row i to row i — Ak (It should be recalled that, 
according to our labelling scheme, we number the rows of M-^ from to 16A;^ ~ 1-) 
Thus, we obtain the determinant of the following matrix: 



fXi X2 X3 X4fe \ 

^1 ^2 

Gi F[ FLk-i 

G2 



G. 



\ 



Fik-2 



F^ 



G^k-i F[ 
• . G^k I 



where the x^'s, I = 1,2,..., 4A;, are defined as earlier, and where the {F'()^s and G^'s, 
i= 1,2,. .. ,4k, are the (4A; — 1) x (4A;) matrices 



/ {-4)**Zi 

' (-4)*'=Z2 



(-4) 





^4fe-2 



(-4)^ 



and 



/ -4X1 -4(^-l)fe-2 
' -4X2 -4{e-l)k-3 



G, 



\ 



-AX^k-i -itk / 



Next we "make" the submatrices Gt 



1,2,..., 4:k, to diagonal matrices, by 



subtracting (j + l)/AXj times column j from column j + 1, j = 1, 2, . . . , 4/c — 1, 
[j + l)/AXj_^k times column j from column j + 1, j = 4A; + 1, 4A; + 2, . . . , S/c — 1, . . . , 
and (j + l)/4Xj_i6fe2+4fe times column j from column j j — 16A;^ — 4/c+ 1, 16A;^ — 
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4/c + 2, . . . , 16A;^ — 1. After these operations we obtain the determinant of the matrix 



/ Xi 



X2 





-^2,2 


G 







G 









2^3 ^Ak \ 

-^3,3 Fik,4k 

F2,Z -^4fe-l,4fe 

-^1,3 Fik-2,Ak 

G ■-. ; 







V 



G 




F2,4k 
Fl^Ak 

G 



J 



where G and the -Fq./j's, 1 < a < P < Ak, are the {Ak — 1) x (4A;) matrices 



/ -4Xi ... 
/ -4X2 



G^ 



-4X4fc_2 . 
-AX^k-iO/ 



and 



^ Jl,2 /l,3 







,(a,/3) ,(a,/3) 
J2,3 J2,4 



•'l,4fe 

f(«,/3) 

•'2,4fc 



Aa.0) „(a,/3) 
.'4fc-2,4fe-l JAk-2,ik 

n f(«./3) 

" •'4fc-l,4fc . 



where 



(",/3) 
rs 



-4) 



(4(/3 - + r + 2), 



Now we ehminate the last columns in F^p for 1 < a < /3 < 4A;. We start by 
eliminating the last column of Fi^^k- We do this by adding f^^l^^ /AXr times column 
16fc^ — 8/c + r to column 16A;^, r = 1,2, ... ,4k — 1. This makes all the entries in 
the last column which are in rows 16A;^ — 8k + 2, . . . , 16A;^ — 4k — 1, 16k^ — 4k zero, 
whereas the entries in the last column in rows 16k^ — 12k + 3, . . . , 16k^ — 8k + 1 are 
modified. Next we eliminate these entries in a similar fashion, by using the columns 
— 12k + r, r — 1,2, . . . ,4k — 1, etc. In the end all the entries in the last column 
in rows 4k, 4k + l, . . . , — 1 will be zero, whereas the entries in the last column in 
rows 1,2, ... ,4k — 1 will have been (significantly) modified. An analogous procedure 
is applied to eliminate the entries in the last columns of Fa,4fc-i, a = 1, 2, . . . , 4A; — 2. 
Just to mention the first step: We add f^^^^'' ^V^^r times column IG/c^ — 12k + r to 
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column IGk"^ — 4k, r = 1,2, ... ,4k — 1. This makes all the entries in column 16k'^—4k 
which are in rows IGk"^ — 12k + 3, . . . , IGk"^ — 8k, IGk"^ — 8A; + 1 zero, whereas the entries 
in the last column in rows 16A;2 -16k + 4,..., IQk^ -12k + 2 are modified. Etc. 

The advantage after having done all this is that now all the entries in columns 
4ik, i = 1,2,..., 4k, are zero except for entries in rows 0,1, ... ,4k — 1. This fact, 
and the fact that the submatrices G are diagonal matrices (of rectangular form) with 
last column consisting entirely of zeroes, makes it possible to reduce the determinant 
of the (new) matrix significantly. For i = IQk"^ — 1, 16/c^ — 2, . . . , 4fc we may expand 
the determinant with respect to row i, in this order. If the details are worked out, 
then we see that our original determinant det is equal to 



4fc-l 



U2 

ei,2 



(_4)i6fc^-4fc -Qxf det 

a=l 

\e4fc-l,l €.4k-l,2 

where is given by ( |4.5|) , and where Cij is given by 



U4k \ 



e4fc-i,4fc/ 



r=l 4A;=to>ti>->tr=i 

j+l=rag>ni>ra2>->nr=l 



-iy-\-4y>^-^''+'+^Zt,Zt,---Zt.^ 



UZlWnu-l)-4k + U+, + 2) 



X 



1,2,...,4A;, and 



Clearly, we may extract (—4)-''^ from column j, j 
(— 4)"^^''"*+-'^/Xj+iXj+2 • ■ • X^k-i from row i, i = 1,2, . . . ,4k — 1 (still using our non- 
standard labelling scheme where the rows are numbered 0, 1, . . . , 4fc — 1), so that we 
obtain the expression 



Ak-l 



(^_4)16fc2-4fe+2A;2(4fc+l)-(2A;-l)(4fc-l) JJ' j^ik+l- 



a=l 



X det 



/ Ul U2 

eii ei 2 



M4fe \ 
GlAk 



\e4fc-i,i e4fe_i,2 ••• e4fc_i,4fc,y 



(4.16) 



where ue is given by 

(_l)^-i(_4)fc8fc(4fc + 1) I n (4^^ - 1) ) I 11(4'^ + 1^ 



^4fc- 



1=1 
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and where Cij is given by 



E E i-iy-'z,,z,,...z,^ 



r=l 4fc=to>ti>->ir=« 

j+l=rio>ni>ra2>->rir=l 



r-1 

X 

iy=0 



l[{Ak{n, -1)-Ak + + 2)t_i,+,_i. (4.17) 



From ( [4.16| ) it is abundantly clear that na=i^ ^t^^^ " divides det (which, after 
all, is equal to ( [4.16| )). It remains to show that also na=7^nfe=o(^a ^ divides 
det M^. 

Let a and 6, l<a<4A; — 1,0<6< a — 1, be given. We want to show that Za—^hk 
divides det M^. We will show the equivalent fact that the rows of det M^\^ are 
linearly dependent. The crucial observation, from which this claim follows easily, is 



that the entries ejj, j = 1, 2, . . . , 4/c, in row i of the determinant in ([4.16|) are given 
by a polynomial in j, Pi{j) say, of degree {Ak — i — 1) and with leading coefficient 
(i.e., coefficient of j^fc-j-i-j gq^g^i 

± Y. (-^r-'Z,Z,--Z., JpT" ■ (4.18) 
This is seen as follows. The summand in (|4.17|) is a polynomial in j, ni,n2, ■ ■ ■ , Ur-i 



of multidegree {Ak — ti — l,ti—t2 — l,t2 — ts — l,.--, U-i —tr — 1) (i.e., the degree in 
j is ik — ti — 1, the degree in hi is ti — 12 — 1, etc.). Because of the fact that (for fixed 
u and varying v) Yl^=u 7^ ^ polynomial in v of degree e + 1 with leading coefficient 
(i.e., coefficient of v'^'^^) equal to l/(e + l), successive summation over hr-i, hr-2, • • • , 
hi yields the claimed facts. 

Under the specialization Za = 4:bk, it is seen by "inspection" that ([4.18|) with 
i replaced by a — 6 vanishes, because the summand corresponding to to > • • ■ > 
ttj-i > C = a > tuj+i > ■ ■ ■ > tr cancels with the summand corresponding to 
to > • ■ • > tuj-i > t^^+i > • • ■ > tr- Hence, the polynomial Pa~b{j) has degree (at 
most) Ak — a + b — 2 (instead of 4A; — a + 6 — 1). Consequently, if Za = Abk then the 
entries in rows a — b, a — 6 + 1, 4k — 1 are given by polynomials in j (to wit: 
j denoting the column index of the entries) of respective degrees 4k — a + b — 2 (!), 
4k — a + b — 2, 4k — a + b — 3, . . . , 2, 1. These are 4k — a + b polynomials, all of degree 
at most 4k — a + b — 2. It follows that there must be a nontrivial linear combination 
of these polynomials that vanishes. Hence, the rows a — b, a — b + 1, . . . , 4k — 1 are 
linearly dependent, which, in turn, implies that the determinant in ([4.16|) (and, thus, 
also det M^) vanishes for Za = 4bk, l<a<4k-l,0<b<a-l. 
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Step 2. Comparison of degrees. Clearly, the degree of det as a polynomial 
in the Xf's and Z^'s is at most IQk^ — 1, whereas the degree of the product on the 
right-hand side of ( [4.15 ) is exactly 16A;^ — 1. Hence, we have 



4fc-l a-1 



det = C2 n [K'^'"" - 4^^)) ' 



a=l b=0 

where C2 is a constant independent of the X^'s and Zts. 

Step 3. Computation of the leading coefficient. In order to determine C2, we deter- 
mme the coefficient of Hali ^^""^a in the expansion of det . By arguments 
similar to those at the end of the proof of Theorem ^, it is seen that this coefficient 
is given by the determinant of the following matrix, which we denote by . It is 
defined exactly in the same way as (see ( [4.9|) ), except that the definitions of the 
functions fo, fi,go, gi are replaced by 

/o(t,j) = (-4f, 

^(t,J) = 4(-4)^ 

9oit,j) = 0, 
gi{t,j) = -A. 

By expanding this determinant with respect to row 0, we obtain 

4fc 

det = ^(-l)^-^M^ det M[, (4.19) 

e=i 

where ue is, as earlier, given by (|4.5| ), and M[ is the matrix arising from by 
deleting row and column 4ik. 

Let i, 1 < £ < 4k, be fixed. We will next compute detM[. When we built M[ 
from M^, we deleted in particular row 0. Therefore we will now switch to the usual 
labelling scheme for rows and columns of a matrix, i.e., we will subsequently not only 
label the columns by 1,2,... but also the rows. 

Ifi < 4k, then we expand det M[ with respect to the last 4k — 1 rows and then with 
respect to the last column. Since these rows and this column contain only one nonzero 
entry, we obtain some multiple of the determinant of a (16/c^ — 4fc — 1) x {lQk^—4k — l) 
matrix. If £ < 4A; — 1, then we continue by expanding the (now reduced) determinant 
with respect to the last 4fc — 2 rows and then with respect to the last 2 columns. We 
continue in the same manner until we have reduced det M[ to the determinant of a 
{4ik — 1) X [4£k — 1) matrix, more precisely, until we arrive at 

detM^^ = 4C'2')"(2)(-4)K''t^') detM;, 
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where M« is the matrix 



(F 












0\ 


G 


F 















G 


F 

















G 


F 















G 


F' 
















{i occurrences of F), 



with F and G the (4/c - 1) x (4 A;) matrices 



4(-4)^^ (-4)^- 



F = 



4(-4)'^ (-4)^^ 

4(-4)'= (-4)V 



and 



/-4 
0-4 



G^ 



0-4 
V 0-4 0/ 

F' the {Ak — 1) x (Ak — 1) matrix which arises from F by deleting its last column, 
and U the (£ - 1) x {4k - 1) matrix 



/-4 





-4 



0\ 




\ -4 . 

We continue by expanding with respect to the last i — 1 rows. Thus we obtain 
det = 4('2')-(2)+^-i(-4)K''t+') det Me (4.20) 



SOME NEW FORMULAS FOR tt 



25 



for det M/', where is the matrix 



(F 








G 


F 








G 


F 








(£ — 1 occurrences of F) , 



... G F 
\0 G VJ 

with V the (4A; - 1) x (4/c - £) matrix 

/ 





(-4)'= 

4(-4)'= (-4)*^ 

4(-4)*^ 



\ 







-4)*^ 



V 



4(-4)'= (-4)^= 
4(-4)V 



(4.21) 



Next we prepare for a reduction from the top of M^. We subtract 4 times column 
j from column j - 1, j = A{i - l)k, 4(£ - l)k - 1, . . . , Ak{i - 2) + 2, j = Ak{i - 
2), Ak{i - 2) - 1, ... , Ak{£ - 3) + 2, j = Ak, Ak - I, . . . ,2, in this order. Thus 
det is converted to det M^", where M^' is the matrix 



/F" 








o\ 


G" 


F" 











G" 


F" 











G" ■•. 





















G" F" 



G" 
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with F" and G" the (4A; - 1) x {Ak) matrices 



/O (-4)^^ 
(-4)'= 



and 



\ 

( -4 

(-4)^ 



_4)4fc-2 




-4 



(-4)*^ 
(-4)V 



;_4)3 



-4 
-4)2 -4 0/ 



Our next goal is to "push" the (nonzero) entries in columns 4tA;-|-l, t = 0,l,...,£— 
2, down to rows (4A; - 1)(£ - 1) + 1, . . . , (4A; - 1)£ - 1, (4A; - 1)^. (This is similar to 
what we did in Step 1 when we "pushed" all the nonzero entries in columns 4i/c, 
t — 1, 2, . . . , 4A; up to rows 0, 1, . . . , 4A; — 1.) In order to achieve this for the 1-st 
column, we add 

J] ^(-l)-(-4)-'-'=+^-i ( - 1 ) ■ (co^^^^ (4^^ + ^)) 
to column 1. Similarly, in order to achieve this for the (4A; + l)-st column, we add 

1-2 4k 



^(_l)^-i(_4)-('-i)fe+s-i 



r=2 s=2 



r-2 



■ (column (4/cr + s)) 



to column Ak + 1. Etc. As a result, the determinant detM^' is converted to the 
determinant of the matrix 



fF" 

G' F" 

G' F" 

G' 












G' F" 

\Hi H2 He_2 G" V) 



(4.22) 
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where F" and G" are as before, G' is the (4 A; — 1) x (4/c) matrix 



G' 









Vo 




-4 



'_4)4fc-3 
'_4)4fc-2 








\ 



-4)^ 
;_4)3 



-4 
-4)2 -4 0/ 



and Hti t = 1, 2, — 2, is a (4A; — 1) x (4A;) matrix with all entries equal to 0, 
except for the entries in the first column. To be precise, the entry in the first column 
and row s of Ht is given by 



-1 



-4)- 



-{l-l-t)k+s 



S - 1 
-l-t 



(4.23) 



It should be noted that the entries in the first column of G" {G" appearing at the 
bottom of the matrix ( |4.22| ), as do the matrices Ht) is given by ( [4.23| ) with t = ^ — 1. 

Now everything is prepared for the reduction. We expand the determinant of 
(|4.22| ) with respect to rows 1, 2, ... , (4A; — 1)(£ — 1). This reduces the determinant of 
(S) to 

(_l)(0(_4)(^-i)(4fc-i)fc detM,, 
where Mn is a (4A; — 1) x (4A; — 1) matrix of the form 

{N V)^ 

with the (s,t)-entry of being given by (|]2D, s = 1, 2, . . . , 4fc- 1, t = 1, 2, . . . , £- 1, 
and V the (4A; — 1) x (4A; — t) matrix from above. If we substitute all this in ( [4.20|) , 
we obtain that det M/' is equal to 



(-l)(D4(")-(^)+^-i(-4)K''t+^)+(^-i)(4fc-i)fc det M,. 



(4.24) 



The submatrix V of Mi is almost diagonal. Subtraction of 4 times row s from row 
s + 1 in Mi, s = i — . . . ,4k — 2, will transform it into a completely diagonal 
matrix (namely into the matrix on the right-hand side of (^4.21 ) with all entries 
4(— 4)'^ replaced by 0). As a side effect, this will turn the (4/c — 1, l)-entry of M^ into 



-4) 



-(^-2)fe+4fc-l 



Ak - 1 

e-1 
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As is easily seen, the determinant of the in this way modified matrix, say, is 
(_l)4fe-^+(^2')(_4)(4fe-^)fe . ((4^ _ i).entry of M;) 

^W{{sA- s)-e^iiy of M;), 

s=l 

or, exphcitly. 

Substitution of the above in ( [4.24|) yields that the determinant det is equal to 



-i\ fAk - 1 
- 1 



- 1 

Now we substitute this in (|4.19| ). We obtain that det is equal to 

^fc-(4A: + l)(4A:-l)!5: - ^ 



'_-|^\fc2l6A:3+20fc2+14fc-l;„4fc|' 



The sum is readily evaluated by means of the Chu-Vandermonde summation (see e.g. 
1^, Sec. 5.1, (5.27)]), so that we obtain 



-2 
4:k - 1^ 



Since the coefficient of 11^=1^ ^a'^^^ ^^^^ expression on the right-hand side of 

( |4.15| ) is exactly the same, we have completed the proof of the lemma. □ 
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